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Abstract 

A method for calculating turbine stage per- 
formance is described. The usefulness of the 
method is demonstrated by comparing measured and 
predicted efficiencies for nine different stages. 
Comparisons are made over a range of turbine pres- 
sure ratios and rotor speeds. A quasi-3D flow 
analysis is used to account for complex passage 
geometries. Boundary layer analyses are done to 
account for losses due to friction. Empirical 
loss models are used to account for incidence, 
secondary flow, disc windage, and clearance 
losses. 
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Nomenclature 

parameter in profile loss calculation 
rotor clearance, fraction of span 
coefficient in disc windage loss model 
coefficient used in Eq. (1) 
coefficient used in Eq. (1) 
chord, m 

specific heat, J/kg*K 
parameter in profile loss calculation 
kinetic energy loss coefficient 
friction factor 

gap between rotor disc and stationary 
surface, m 
blade height, m 

incidence (gas minus blade) angle 
kinetic energy loss, J/kg 

flow rate, kg/sec 
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pressure, N/m 
exponent used in Eq. 
radius, m 
Reynolds number 
rotor reaction 
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loss coefficient 
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blade loading coefficient 
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inlet air angle 
exit air angle 
mean air angle 
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COS a„ 


+ tan a. 


“Z / 

inlet blade angle 

specific heat ratio 

displacement thickness 

blade trailing edge thickness, m 

efficiency 

momentum thickness, m 
viscosity, m-sec/kg 
density, kg/rtW 
kinetic energy thickness, m 
pitchwise angle 

angular distance for one pitch 
speed of rotation, rad/sec 


Subscripts 

oa overall 

cL clearance 

cr critical velocity 
d rotor disc 

e endwall 

ex exhaust duct 

fs free stream 

h hub 

i incidence 

m mixed out location 


■'r 

S 

pitchwise distance between blades, m 

0 

total (ab; 

T 

temperature, K 


rotors) 

V 

absolute’ velocity, m/sec 

’ p 

prof i 1 e 

W 

relative velocity for rotor, absolute 

rot 

rotor 


velocity for stator, m/sec 

s 

secondary 

X 

axial distance between rotor and measuring 

St 

stator 


station, m 

. t 

tip 



w 

windage 
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0 zero clearance 

1 inlet to blade row 

2 exit of blade row 


Superscripts 

* normalized by pitchwise distance 

' absolute total 

" relative total 

Introduction 

In the effort to improve turbine perform- 
ance greater use is. being made of complex blade 
passage geometries. I 1- '*' These geometries can 
include contoured endwalls, and blades that are 
leaned and bowed, resulting in nonuniform work 
distributions. A method for predicting turbine 
stage performance has been developed which 
accounts for the effects of complex passage geom- 
etries. The method uses a quasi-3D inviscid flow 
analysis iteratively coupled to calculated losses 
so that changes in losses result in changes in the 
flow distribution. In this manner the effects of 
both the geometry on the flow distribution and the 
flow distribution on losses are accounted for in 
the prediction of stage performance. This report 
consists of a description of the method, and a 
comparison between predicted and measured stage 
performance for nine different single stage axial 
turbines. The turbines had diameters between 12 
and 75 cm, and design pressure ratios between 1.8 
and 3.1. The present method was developed pri- 
marily as an aid in understanding why a particular 
stage performs as it does. The reason for under- 
taking this effort is that existing correla- 
tions do no t account for all geometric 

factors which affect stage performance. This work 
is an extension of the work described in Ref. 7 
which was used to predict the change in stator 
performance due to endwall contouring. 

The analysis uses a quasi-3D flow analysis 
coupled to a boundary layer analysis to determine 
losses due to friction on the blade surfaces and 
passage endwalls. The flow analysis assumes in- 
viscid flow, and uses existing programs. '°>9) 

The boundary layer analysis also uses an existing 
program.!^' The profile losses for both the 
stator and rotor are determined by using the 
results of the boundary layer analysis along with 
flow and geometry information in the loss model 
of Ref. 11. Empirical models are used to account 
for incidence, disc cavity, rotor tip clearance, 
and secondary vortex losses. In order to demon- 
strate the validity of the method, comparisons are 
made between predicted and measured stage effi- 
ciencies over a range of pressure ratios and rotor 
speeds for each of the nine turbines. A range of 
operating conditions was examined for each turbine 
in order to assess the adequacy of the individual 
loss models. For example, decreasing speed at 
constant pressure ratio increases rotor incidence 
loss. The ability of the analysis to predict the 
change in efficiency when speed is reduced is a 
measure of the adequacy of the incidence loss 
model. The adequacy of the tip clearance model is 
assessed by comparing predicted and experimental 
stage efficiencies at different rotor clearances 
for a single turbine. In addition to stage effi- 
ciencies the analysis provides other data of 


interest to a designer. Comparisons are made 
between predicted and experimental flow angles 
for three turbines at their design conditions. 


Method of Analysis 

The predicted stage performance was calcu- 
lated by first determining the stator performance. 
These results determined the flow field entering 
the rotor. After the rotor performance was deter- 
mined, the overall stage performance was then 
calculated. For each blade row the flow field was 
determined using a quasi-3D flow analysis. This 
analysis was done using existing stream function 
flow analysis codes. First a hub-to-shroud solu- 
tion was obtained using the program MERIDL de- 
scribed in Ref. 8. One of the results of the 
MERIDL calculation is the variation of stream- 
sheet thickness with streamwise location for 
various spanwise locations. The stream-sheet 
thickness is the thickness necessary to pass a 
fixed fraction (one percent in this program) of 
the total flow. These results and other data 
were used as input to obtain the bl ade-to-bl ade 
solution. This solution was obtained using the 
program TSONIC described in Ref. 9. The flow 
conditions at the exit of the blade row were not 
imposed as part of the input. They were iterated 
until the static pressure near the trailing, edge 
was the same for both the pressure and suction 
surfaces of the blade. 

The flow analysis was done using modified 
versions of the MERIDL and TSONIC programs. Modi- 
fications to the program MERIDL are described in 
Appendix A, and TSONIC contained the modifications 
of Ref. 12. These modifications allowed solutions 
to be obtained in the transonic flow regime. How- 
ever, the presence of .shocks was not accounted 
for. Appendix A also describes two additional 
modifications made to TSONIC. One allows blade- 
to-blade solutions to be obtained assuming trail- 
ing edge ejection of sufficient flow to give no 
velocity defect behind the trailing edge. The 
second modification allows the specification of ■ 
an elliptical blade leading edge in the program 
TSONIC. These two modifications are presented to 
document the current status of the program, but 
were not used in the subject investigation. 

Once the flow analysis was done for the 
stator row, boundary layer analyses were per- 
formed. These analyses were done along streamline 
paths for both the blade surfaces and endwalls. 

The analyses were done using the program BLAYER 
described in Ref. 10. Boundary layer parameters 
were obtained along streamline paths determined 
from the inviscid analysis. The profile loss was 
calculated using the method of Ref. 11, and the 
endwall loss using the method of Ref. 7. Appendix 
B gives a summary of the equations used to deter- 
mine losses, and shows how these values were used 
to determine stage efficiency. The loss in total 
pressure used in the flow analyses could be a 
function of span. To do so requires assumptions 
regarding the spanwise distribution of loss. The 
simplest assumption was that there was no spanwise 
variation, and this assumption was used for most 
of the comparisons. After the losses for the 
stator were determined they were compared with the 
assume'd values used in' the flow analysis. If 
there was a difference new losses were assumed for 
use in the flow analysis and the analysis was 
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repeated. This is shown in Fig. 1 which is a 
block diagram of the procedure used to determine 
stage efficiency. 

Once the losses for the stator row had been 
determined, the process was repeated for the rotor 
using the results of the stator analysis as input 
to the rotor analysis. The procedure for calcu- 
lating losses was the same for the rotor as for 
the stator, with the losses being calculated in 
the relative frame. Only hub friction losses were 
calculated for the rotor, since the tip friction 
loss would be included as part of the tip clear- 
ance loss. The duct loss is the loss due to fric- 
tion on the endwalls between the rotor trailing 
edge and the measuring station. A friction coef- 
ficient was used as a simplified boundary layer 
analysis to determine this loss. This was done 
primarily as a convenience, since the program 
BLAYER could have been used to determine this 
1 oss. 


The losses due to incidence, windage, and 
clearance are summarized in Ref. 13. The incid- 
ence loss is based on the assumption that the 
kinetic energy in the normal component of velocity 
is lost. The actual equations used for incidence 
loss incorporate the empirical results of Ref. 14. 
The results of Ref. 14 give a six degree bias in 
the zero loss incidence angle, and a greater loss 
for positive incidence. The disc windage loss 
accounted for the friction between the rotating 
disc and stationary surfaces. The equations of 
Ref. 15 were used. This loss was not applied to 
those stages where the disc windage loss was 
included as part of the tare loss. The loss 
resulting from the clearance between the rotor 
tip and the casing used the data of Ref. 16. The 
change in stage efficiency due to clearance is a 
function of the rotor gap and rotor reaction. 

The equations used to express this functional 
relationship are the same as those given in 
Ref. 17. 

The calculations were carried out on an IBM 
370 computer. Approximately 3 min of CPU time 
were needed to converge a case for each blade row. 

Secondary loss model . The correlation for 
the loss due to the secondary vortex was the only 
loss correlation which was modified for the pre- 
sent analysis. There are two widely quoted sec- 
ondary loss correlations.. These are Morris and 
Hoare/ 1 ®' and Dunham. The equations are 
of a similar form and both can be expressed as: 

h (/) ' * <*] d) 


_ cos “ 2 (c\ 
s cos 8| \hy 


Neither correlation accounts for variations in 
this parameter, while the endwall boundary layer 
analysis did. Therefore, the secondary loss model 
used herein was only the first term of Eq. (1). 

The question of the appropriate coefficients for 
this term remained to be resolved. Figure 2 shows 
a comparison of the two correlations, along with 
the experimental data that was used for the 
correlation of Dunham. The total loss values for 
both correlations are in good agreement. However, 
they differ significantly for the loss due only to 
the inlet boundary layer. The equation of Morris 
and Hoare gives nearly a lower bound for the loss 
due to the inlet boundary layer. Also shown in 
Fig. 2 is a curve for the loss due to the inlet 
boundary layer that is 1.5 times that of Morris 
and Hoare. It is in better agreement with the 
data for those cases where the experimental loss 
either did not or may not have included the end- 
wall loss. Coefficients of Ci = 0.44, C 2 = 0, 
and p = 1.0 were chosen for the secondary loss. 
The equation used for the secondary loss due to 
the inlet boundary layer is 1.5 times the value 
used by Morris and Hoare, and is: 


cos 


cos 
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The secondary loss model was then used to 
compare predicted and measured kinetic energy loss 
coefficients for three stators. Table I gives 
this comparison for experimental data from Ref. 

20. These stators were later used in stage tests, 
and the experimental data for the stage tests will 
be discussed subsequently. The data are for 
stators 13 cm in diameter, which have cylindrical, 
conical, and S-shaped contoured tip endwalls. 

This table shows a breakdown of the analytic 
losses for the three stators. The secondary loss 
due just to the inlet boundary layer is approxi- 
mately one quarter of the total loss for all three 
stators. The ratio of the total wall displacement 
thicknesses to the blade chord at the blade lead- 
ing edge varied between 0.008 and 0.013. It can 
be seen from Fig. 2 that this is in a region where 
the correlation of Dunham and the correlation of 
Morris and Hoare for the secondary loss due to the 
inlet boundary layer significantly disagree. .The 
agreement between the analysis and the measured 
data is good because the secondary loss was chosen 
as 1.5 times the value from the correlation of 
Morris and Hoare. Also, a preliminary screening 
of the data for the nine turbines showed that the 
unmodified equation of Morris and Hoare led to a 
higher than measured efficiency, while the 
equation of Dunham led to a lower than measured 
efficiency. 


Comparisons With Experimental Data 


For the correlation of Morris and Hoare = 
0.294, Co = 0.011, and p = 1.0. For Dunhams's 
correlation = 0.078, Cg = 0.055, and p = 

0.5. The secondary loss is the sum of two terms. 
The first term is a function of the inlet boundary 
layer height, and represents the loss resulting 
from the formation of the leading edge vortex. 

The second term is loss due to boundary layer 
growth on the endwall. This is the same loss as 
is accounted for by endwall boundary layer growth 
calculated by the BLAYER analysis. The endwall 
loss is strongly dependent on the distance between 
the trailing edge and the measuring station. 


Stage descriptions . Some characteristics of 
the nine axial flow turbines are shown in Table 
II. The individual stages are labeled according 
to their tip diameter. A wide range of turbine 
sizes is represented. With one exception all of 
the stages were tested with an inlet temperature 
near ambient, and that one had an inlet tempera- 
ture of only 422 K. The tests were run with 
either the inlet or exit pressure near atmos- 
pheric. The stages are ordered according to their 
equivalent, flow.. Some stages had common blade 
rows. Stage 42A and 42B had the same stator, but 
different rotors. This was also true of stages 
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25U and 25T. These last two turbines differed in 
that the untwisted rotor of 25U was twisted for 
turbine 25T. Stages 13-Cyl and 13-Cone had dif- 
ferent stators, but the same rotor. Turbine 
13-Cyl had cylindrical stator endwalls, while 
13-Cone had a conical tip endwall. 

Stage efficiency comparisons . Figures 3 to 
11 give comparisons of the predicted and measured 
stage efficiency for the nine stages. In each 
figure data are given for three speeds, between 
80 and 110%of equivalent design speed. For each 
stage the analysis was able to predict the effi- 
ciency within 1.2 points at design pressure ratio 
and speed. Off-design comparisons were made for 
each of the nine turbines at constant work factor, 
which is near the operating line. These compari- 
sons are shown in table III. Some turbines do not 
have comparisons at off-design speeds, because 
there was no experimental data available at the 
required pressure ratio. With the exception of 
turbine 51 at 80%speed, the analysis was able to 
predict the efficiencies within 2 points. 

Table IV gives comparisons between measured 
and predicted efficiencies at each design pressure 
ratio for each speed. The analysis was able to 
predict the efficiencies at off-design speeds 
about as well as it did at 100% speed. At 100% 
speed the analysis overpredicted the stage effi- 
ciency by an average of 0.2 point. At 110%speed 
the average overprediction was also 0.2 point. 

At 90% speed the average overprediction was 0.4 
point, but there were only three comparisons. At 
80% speed the average overprediction was 0.1 
point. When the rotor speed is varied, but the 
pressure ratio is held constant, there are signif- 
icant changes in the loss due to incidence. Other 
losses also change, but to a lesser extent. The 
results of the comparisons in table IV indicate 
that the incidence loss model is reasonably 
accurate. 

Figures 3 to 11 show that as the pressure 
ratio is decreased at constant rotor speed there 
was a greater discrepancy between the analysis and 
the experimental data. This is a result of the 
procedure used for the blade boundary layer cal- 
culations. When separation was indicated, the 
surface velocities were smoothed until the cal- 
culation went as far as the trailing edge. For 
many of the rotors, pressure surface separation 
was indicated at low pressure ratios. For 
example, at pressure ratios less than 2.2 the 
analysis showed separated boundary layers at all 
speeds for turbine 51. The analysis showed sepa- 
rated boundary layers, but could not predict the 
penalty associated with the separation. If the 
penalty could be predicted, the efficiency would 
be lower than the efficiency predicted with a 
smoothed velocity profile. For each of the nine 
turbines predictions are shown for pressure ratios 
greater than the design pressure ratio. The 
agreement with the experimental data is about as 
good as it was at the design pressure ratio. The 
predictions do not cover the entire pressure range 
for all turbines. The analytic predictions were 
made only in the region where the rotor exit Mach 
number was less than 1.3. Shock loss would be 
significant only at higher Mach numbers. Turbines 
25U and 25T show a rapid decrease in experimental 
efficiency at the highest pressure ratios. This 
may have been the result of shock losses. 


Clearance effects . Figure 12 shows stage 
efficiency versus rotor clearance for turbine 
13-Cyl at design speed. Results are shown at both 
constant work and constant pressure ratio. In 
each case two predictions are shown. One predic- 
tion assumes rotor loss is uniformly distributed 
from hub to tip. The other prediction has a 
linear loss distribution. In this second predic- 
tion losses other than tip clearance loss are 
uniformly distributed. The tip clearance loss 
varies linearly, with none of the loss at the hub 
and the average loss at mid-span. The turbine had 
a relatively large design clearance of 2.4% of 
span. At the design clearance both loss distribu- 
tion assumptions result in the same efficiency. 

The difference between the predicted efficiencies 
for the two loss distribution assumptions in- 
creases with increased clearance. Even the uni- 
form loss distribution model gives a satisfactory 
prediction for the change in efficiency with 
clearance. 

Loss breakdown . Figure 13 shows the break- 
down of losses for the different turbines at the 
design point. The individual losses are normal- 
ized by the overall loss, so that they are a frac- 
tion of the total loss. The effect of turbine 
size can be seen from the data in this figure. 

The 25 cm turbines and the 13 cm turbines both had 
clearances of nearly 2.5%of span. For the 25 cm 
turbines the clearance loss is about half the 
total loss for the stage, while for the 13 cm tur- 
bines it is only about a third of the total loss. 
Also evident is the growth of the importance of 
rotor secondary flow loss as the size decreased. 
The rotor secondary loss is large because of the 
large amount of flow turning and the thickness of 
the inlet boundary layer. In reference to Fig. 2, 
the total inlet boundary layer displacement 
thickness for turbine 13-Cyl was only 0.01 of the 
rotor chord. This turbine has the largest loss 
due to rotor secondary flow of the nine turbines. 
At this value there is a large difference between 
the correlation of Morris and Hoare and the cor- 
relation of Dunham for the secondary loss due to 
the presence of the endwall boundary layer at the 
blade leading edge. The rotor secondary loss is 
an uncertain loss. Added to the uncertainty 
arising from the cascade data is the uncertainty 
arising from skewing of the rotor inlet boundary 
layer. A 50% variation in the rotor secondary 
vortex loss would account for most of the differ- 
ence between the predicted and experimental effi- 
ciencies. 

Flow angles . A comparison between measured 
rotor exit flow angles and flow angles predicted 
by the analysis is given in Fig. 14. The compar- 
ison are for turbines 12, 13-Cyl, and 13-Cone at 
the design condition. Viscous effects play a 
significant role in determining the local flow 
angle. The predictions are reasonably good con- 
sidering that the analysis has no spanwise varia- 
tion in loss, and that the flow angles were deter- 
mined from the quasi-3D inviscid flow analysis. 

Concluding Remarks 

A method for predicting turbine stage per- 
formance has been developed. This method uses a 
quasi-30 inviscid flow analysis iteratively 
coupled to calculated losses to predict stage 
performance. At the design pressure ratio and 
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speed it was able to predict the stage efficiency 
within 1.2 points. The agreement was not as good 
at other conditions. Near the operating line the 
analysis was generally able to predict the effi- 
ciency within 2.0 points. At low pressure ratios 
and design speed the analysis generally predicted 
a significantly higher efficiency than was meas- 
ured. The analysis indicated when separation was 
likely to occur, but proceeded to smooth the sur- 
face velocities until no separation occurred. 

This was done before the loss prediction was made. 
For the cases where separation occurred, it could 
be inferred that the actual efficiency would be 
less than predicted. Because of the high losses 
associated with separated flow it is desirable to 
be better able to predict the condition at which 
separation begins, and to be able to calculate the 
losses when separation occurs. To improve the 
precision in determining when separation occurs, 
if it occurs near the leading edge, a blade-to- 
blade solver such as described in Ref. 26 would be 
useful. This program determines a more accurate 
solution in the leading edge region through the 
use of body fitted coordinates. In order to 
better define profile loss at other than design 
conditions, a method of calculating loss with a 
separated boundary layer would be desirable. How- 
ever, this was beyond the scope of this work. 

The analysis predicted stage efficiencies at 
high pressure ratios reasonably well. The analy- 
sis was not used in regions where shock losses 
would be significant. Use of flow analyses which 
would account for these losses, as well as 
directly accounting for viscous losses, would 
increase confidence in the model. However, use of 
codes which incorporated these features would have 
resulted in much larger computational times. This 
in turn would have lessened the number of cases 
for which comparisons could be made. 


Appendix A - Description of Modifications 
to the Programs MERDIL and TSONIC 


Trailing Edge Ejection . One of the diffi- 
culties in using the TSONIC program^) is 
determining the outlet flow angle, BETAO. (This 
appendix uses the nomenclature of Refs. 8, 8a, 
and 9.) If the trailing edge was sharp, there 
would be no difficulty, since the Kutta condition 
could be used. However, on real turbomachines 
blades, there is considerable thickness at the 
trailing edge. The usual practice has been to 
visually examine the velocities near the trailing 
edge on the suction and pressure surfaces, and 
then to adjust BETAO until the velocities on the 
two surfaces approach a common value at the trail- 
ing edge. This usually requires good judgement 
because of the velocity peaks induced by the flow 
going around the small trailing edge radius. An 
alternative scheme has been devised, whereby the 
trailing edge radius may be replaced by a blunt 
trailing edge. Because of the large blockage, and 
to avoid flow around a sharp corner, flow is con- 
sidered to be issuing from the trailing edge at 
the average velocity of the freestream flow 
between the blades at the trailing edge. This 
also has the effect of compensating for wake 
blockage downstream of the blades. 


The blade configuration at the trailing-edge 
should have a small uniform curvature near the 
trailing edge with a moderate wedge angle (see 


Fig. 15). If the wedge angle is small the veloc- 
ities will vary smoothly near the trailing edge. 
With a larger wedge angle there will tend to be a 
reduced surface velocity at the trailing edge. 

In either case the trailing edge velocities will 
vary smoothly with a variation of BETAO. This 
makes its feasible to write coding to vary either 
the downstream whirl or BETAO to obtain closure of 
the suction and pressure surface velocities. 

Four subroutines are affected by the changes 
required to simulate trailing edge flow ejection. 
First, the blade geometry is modified to have a 
blunt trailing edge. This is done in BLCD. 

Then, because of the additional mass flow from 
the trailing edge, the downstream periodic 
boundary condition is changed to allow for the 
increased mass flow. This requires changes to 
subroutines COSUB and VELSUB. Another change 
that is required is to provide for proper bound- 
ary values along the blunt trailing edge (see 
Fig. 15). This requires changes to COSUB and 
VELTAN. 


Figure 15 shows the geometry of the blunt 
trailing edge as derived from the original trail- 
ing edge radius. The tangent line to the trail- 
ing edge radius is extended to the trailing edge 
as shown in Fig. 15. The equation for the dis- 
tance from the center of the trailing edge to the 
intersection point is 

The coordinates of this intersection point is used 
for the last spline point on the blade, and 
replaces the spline point at the tangent point on 
the trailing edge radius. 

Figure 15 also shows the mesh geometry at the 
trailing edge. Additional flow comes out the 
trailing edge. It is assumed that the ejection 
velocity is the average velocity across the pas- 
sage at the trailing edge. The stream function 
(u) varies from zero to 1 from the lower blade to 
the upper blade. Then the stream function in- 
creases linearly across the blunt trailing edge 
of the upper blade. The total increase is equal 
to the tangential blade thickness divided by the 
tangential passage width at the trailing edge. 


Elliptical Leading or Trailing Edge s. The 
TSONIC program uses subroutine BLCD to calculate 
the blade geometry. The version of BLCD published 
in Ref. 9 uses circular arcs at the leading and 
trailing edges. This is satisfactory for many 
blades, however, in some cases it is not. One 
way to obtain a more general leading edge shape 
is to use an ellipse which may be at an angle. 

The only subroutine changed is BLCD. The input 
format was not changed, but the definitions of 
several input variables were changed. The new 
definitions are given in Table V, and are illus- 
trated in Fig. 16. See Ref. 9 for the input for- 
mat. The equation of an ellipse at an angle, 
BETI1, to the x-axis is: 


Cx ± RI1 * RI2 
A 


where 


A = RIl 2 cos 2 (BETIl) + R 1 2 2 s i n 2 (BETIl 1 
C = s i n ( BET II) cos (BETI1) (RI1 2 -RI2 2 ) 
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The maximum and minimum value of x is ± ■\JK. 
The corresponding value of y is Cx/A. 
Therefore: 

x • = at 

min V 


x max = V*’ at 

Equations (A2) and (A3) are used in BLCD to cal- 
culate the e coordinate of the blade surface 
for any given m-coordinate, using y + C I -JK = re, 
and x + = m on the el 1 ipse. 

Mesh Generation for MERIDL . The MERIDL 

program^ generates an orthogonal mesh, which 
is used to obtain a finite difference solution to 
the stream function equation. The scheme for 
generating the solution mesh is described in Ref. 
8a. After generating the streamwise mesh lines, 
the hub-to-shroud mesh lines are generated to be 
orthogonal to the streamwise mesh lines. The hub 
is divided into equal increments within the blade, 
and upstream and downstream of the blade. In 
some cases with high wall curvature, this can 
result in crowding of the hub-to-shroud mesh lines 
in some region, and wide distances between lines 
in another region. This can be alleviated some- 
what by starting at a streamwise mesh line near a 
region of high curvature. Figure 17 shows the 
grids generated using different starting loca- 
tions. When the grid is uniform at mid-span, the 
grid is widely spaced in the region of high cur- 
vature. When the grid is started two lines away 
from the curved wall, the spacing is better at 
this wall, and not too compressed at the opposite 
wall. The only subroutine that is affected is 
MESHO. The revised subroutine starts at the mid- 
dle of the passage and goes both directions. The 
index of the initial streamwise mesh line is JMID. 
JMID can be set to any value from 1 to MHTP1, (the 
mesh line at the tip), to determine the starting 
line for generating the hub-to-shroud mesh lines. 

Alternate Method for Calculating Tangential 
Velocity Components . The basic MERIDL program 
uses an empirical correction to the blade shape 
near the leading and trailing edges. This is 
done to avoid a discontinuity in tangential momen- 
tum when there is nonzero incidence or deviation 
angles. The empirical correction is described 
Appendix F of Ref. 8. The correction, however, is 
not satisfactory in many cases. It is particu- 
larly poor in cases with low solidity (1/2 or 
less). In some cases with low solidity, the flow 
does not become parallel to the blade surface 
anywhere in the passage. An alternate scheme was 
devised to improve the calculations for cases with 
very low solidity. Instead of using the blade 
geometry angles directly, the tangential momentum 
is used. The tangential momentum upstream and 
downstream of the blade is specified as a function 
of the stream function. For the new scheme the 
tangential momentum within the blade passage is 
assumed to vary in proportion to the variation of 
the tangent of the streamwise blade angle. If the 
passage height varies by less than 50% between the 
leading and trailing edges, the effect on the BESP 
array used by TSONIC is not large. 



There is an additional advantage to specify- 
ing the tangential momentum instead of the blade 
angle within the blade. When the blade angle is 
specified, the stream function equation becomes 
hyperbolic when the relative velocity becomes 
supersonic on the mean stream surface. If the 
velocity becomes supersonic anywhere within the 
blade on the stream surface, MERIDL cannot cal- 
culate the solution. However, when the tangential 
momentum is specified, the stream function equa- 
tion does not become hyperbolic until the meridi- 
onal component of the relative velocity becomes 
supersonic. Normally, the meridional component of 
the velocity will be subsonic, even with super- 
sonic relative velocity, so that a solution can 
be obtained. 

Two subroutines were changed and a new sub- 
routine (TANVEL) was added to calculate the 
tangential velocity components. The velocity 
components are calculated by NEWRHO. In the 
original version, W 6 is calculated from the 
meridional components W s and W(- by using the 
blade slopes, de/ds and de/dt. The alternate 
method calculates W 6 by varying the whirl 
(rV e ) in the same manner as the variation in 
rde/ds between the leading and trailing edges. 

This calculation is done by TANVEL, which is 
called by NEWRHO. The rest of NEWRHO is essen- 
tially unchanged. In addition INDEV was changed 
to eliminate the calculation of de/ds near the 
leading and trailing edges. 


Appendix B - Summary of Equations 
Used In Loss Model 


This appendix describes the equations used to 
calculate losses which are in turn used to calcu- 
late stage efficiency. These equations are spe- 
cifically for an axial machine. However, the 
loss calculation is implemented in a general way, 
and the performance of a radial stage can be cal- 
culated using the method described in Ref. 27. 

Profile Loss . - The profile loss is calcu- 
lated according to the method of Stewart. 

It includes the effects of boundary layer growth 
along the blade surfaces as well as the thickness 
of the blade trailing edge (station 1). The pro- 
file loss is calculated at a location where com- 
plete mixing is assumed to have occurred (station 
m). The velocity and total pressure are uniform 
in the pitchwise direction. To calculate the 
kinetic loss coefficient two parameters (C and D) 
are first calculated: 
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where s and 6^ p are the displacement 
thickness’on the suction and pressure surfaces 
near the trailing edge. Values are obtained from 
the boundary layer analysis which terminated one 
grid line upstream of the trailing edge. 

Simi larly. 
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where 6t e is the blade thickness normal to the 
chamber line. The ratio of the component of ve- 
locity in the streamwise direction to the critical 
velocity ( W x / W c r ) m is found from: 



The density ratio (o/p 0 ) m is calculated from: 



The total pressure ratio is calculated from 
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Secondary Loss - The secondary loss is cal- 
culated for the inlet boundary layer on each end- 
wall using equation 1., This equation does not 
account for any variation in free stream proper- 
ties across the span. To account for the varia- 
tion in an annular cascade a mass average of the 
secondary loss is done using the inviscid flow 
properties between the hub and tip. This results 
in: 


[ ( Vfs“fs cos ar) h * (Y s p fs W fs cos ar) t] 


J r h 
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r Pf s Wfs cos a dr 


From the definition of Y the equivalent e can be 
calculated. 


Incidence Loss - The kinetic Toss due to 
incidence is given in Ref. 14 as: 
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The profile kinetic energy loss coefficient is 
calculated from: 
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The mass average profile loss for the entire 
passage is: 


P = 



pW xV dr 


pW^r dr 


(Bl) 


where i op t = -6° and n = 3 if i > i op t and’ 
n = 2 if i < i 0 _ t . Ibis co nver t e d into an 
equivalent temperature drop, which in turn is con- 
verted into an which represents a loss in 
total pressure due to incidence. 



Tip clearance Loss - The correlation for the 
tip clearance loss requires that the stage effi- 
ciency at zero clearance (no) be known. The ratio 
of stage efficiency to the efficiency at zero 
clearance is: 

= 1 - C, (1.72 + 0.108R + 2.755 R 2 ) 

n 0 r r 

C-| is the clearance fraction and R r is defined 
as: 



Endwall Loss - The endwall loss is determined 
from boundary layer calculations. The boundary 
layer growth is calculated for a number of invis- 
cid streamlines in the blade passage. For the 
stator the streamlines extend from the upstream 
boundary to the rotor leading edge along both the 
hub at tip. For the rotor only boundary layer 
growth along the hub between the leading and 
trailing edges is considered. The endwall loss, 
derived in Ref. 7, is: 


Exhaust Duct Loss - The loss accounts for the 
pressure drop between the rotor exit and the 
measuring section 


4f 


p V 2 
ex ex 
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ex 2 ( r t - r h ) cosa 2 ’ 

friction factor for a smooth tube at the exhaust 
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duct Reynolds number, the flow path length is 
x/cosa, and the hydraulic diameter is 2 (rr - 
r h) • 


Disc Windage Loss - The rotor disc windage 
loss is given is Ref. 15 as: 

2 

»r d 

(Re)d = — , where u = rotor angular velocity, 

r d = the disc radius, v is the dynamic viscosity, 
and p d is the average density at the hub. The 
loss coefficient is a function of the gap (g) and 
the Reynolds number. 

C w = 0.102 0-1 (Re )' 0 ' 2 if (Re) d > lxlO 5 

/ \ -0.167 A «r r 

C w = 0.080 (f-j (Re)-°' 2b if (Re) d < lxlO b 

The disc friction power loss is: 

„ 3 5 

C w 0hl r d 
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The pressure ratios can be expressed in terms of 
stator and rotor kinetic energy loss coefficients. 



(B7) 


The equivalent kinetic energy loss coefficient is: 



Stage Efficiency - The stator kinetic energy 

loss coefficient (e s t) is taken as the sum of the 
loss coefficients for the stator. The rotor 
kinetic energy loss coefficient at zero clearance 

(e ro t o) 1S taken as the sum of the rotor kinetic 
loss coefficients for the incidence loss, profile 
loss, endwal 1 loss, secondary loss, and disc windage 
loss. The stage efficiency is calculated in three 
steps. First the efficiency at zero clearance be- 
tween the stator inlet and rotor exit is determined. 
Second, the change in efficiency due to rotor clear- 
ance is determined using equation B5. Third the 
overall stage efficiency is calculated by including 
the pressure drop in the exhaust duct. Figure 18 
is a sketch of the paths followed in calculating 
the stage efficiencies. P} ' - P 2 1 is the loss 
in total pressure in the stator. P 2 " - P 3 " is 
the loss in total pressure in the rotor excluding 
the clearance loss, while P 2 " - P 4 11 includes 
the clearance loss. P 4 ' - P 5 1 is the loss in 
total pressure in the exhaust duct. The zero 
clearance efficiency is: 



The efficiency including clearance is calculated 
from equation B5 to give: 

n cL = 10 - An c i_ 

Since there is no additional work, but there is a 
pressure drop in the duct. 
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TABLE I. - COMPARISON OF MEASURED AND 
PREDICTED STATOR KINETIC ENERGY 
LOSS COEFFICIENTS FOR THREE 
13 cm DIAMETER STATORS 



Kinetic energy loss 
coefficients 

< s ,t> 

Stator 

13-Cyl 

13-Cone 

13-S shaped 

Predicted 



■■■ 

P rof i 1 e 

0.027 

0.024 

■ 

Endwal 1 

' .017 

.016 

irill 

Secondary 

.010 

.011 

.010 

Incidence 

.001 

.001 

.001 

Total 

.055 

.052 

.052 

Measured^) 

0.055 

0.052 

0.052 






TABLE II. - TURBINE CHARACTERISTICS AT THE DESIGN CONDITION 







Turbine 





Label 

76 cm 

51 cm 

42-U 

42-T 

25-U 

25-T 

12 cm 

13-Cyl 

13-Cone 

Tip diameter, cm. 

76 

51 

42 

42 

25 

25 

12 

13 

13 

Hub/Shroud radius ratio 

Mean radius Reynolds no. (xlO'^' 

0.75 

0.86 

0.85 

0.85 

0.85 

0.85 

0.73 

0.83 

0.81 

2.7 

8.1 

8.4 

8.4 

8.7 

8.7 

2.4 

3.1 

3.1 

Inlet total pressure, atm 

1.0 

3.1 

1.0 

1.0 

1.7 

1.7 

0.8 

1.4 

1.4 

Inlet total temperature, K 

300 

422 

294 

294 

306 

306 

326 

323 

323 

Stator 










Axial aspect ratio 

2.4 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

0.7 

0.7 

Axial solidity , 

Pressure ratio - Pi / ?2 

1.0 

0.6 

0.6 

0.6 

1.1 

1.1 

1.6 

2.0 

1.9 

1.5 

1.9 

1.9 

1.9 

1.6 

1.6 

1.6 

2.0 

2.0 

Rotor 










Axial aspect ratio 

1.9 

1.1 

1.0 

1.0 

1.1 

1.1 

1.2 

0.5 

0.5 

Axial solidity 

1.6 

1.1 

1.4 

1.4 

1.5 

1.5 

2.0 

1.7 

1.7 

Rotor speed, rpm 

4650 

9048 

10 686 

10 686 

12 750 

12 750 

29 440 

33 300 

33 300 

Clearance, % span 

0.75 

0.80 

0.86 

0.86 

2.30 

2.30 

1.20 

2.40 

2.40 

Stage 










Total -to-total pressure ratio 

1.75 

2.72 

3.1 

3.1 

1.82 

1.82 

2.05 

2.77 

2.77 

Reference for experimental data 

21 

22 

3 

3 

23 

23 

24 

25 

25 



















TABLE III. - COMPARISON BETWEEN MEASURED AND PREDICTED EFFICIENCIES AT CONSTANT WORK FACTOR 


Turbine 


% equivalent speed 
Total -to-total pressure ratio 
Measured efficiency 
Predicted efficiency 

% equivalent speed 
Total -to-total pressure ratio 
Measured efficiency 
Predicted efficiency 

% equivalent speed 
Total -to-total pressure ratio 
Measured efficiency 
Predicted efficiency 


76 

51 

42 A 

42 B 

100 

100 

100 

100 

1.75 

2.36 

3.09 

3.09 

.922 

.895 

.920 

.911 

.923 

.907 

.916 

.915 

90 

1.57 

110 

2.91 

no 

no 

.914 

.913 

.895 

.911 

N. A 

N. A 

80 

80 

90 

90 

1.42 

1.71 

2.45 

2.41 

.919 

.857 

.894 

.900 

.917 

.888 

.905 

.908 


25 U 

25 T 

12 

13-Cyl 

100 

100 

100 

100 

1.82 

1.82 

2.05 

2.77 

.874 

.882 

.863 

.845 

.881 

.876 

.865 

.850 

no 

2.10 

no 

2.10 

no 

no 

.885 

.879 

.894 

.881 

N. A 

N. A 

80 

80 

80 

80 

1.45 

1.45 

1.57 


.865 

.890 

.857 

N. A 

.874 

.883 

.851 



13-Cone 


100 

2.77 

.855 

.853 

110 


N. A 


80 


N. A 


N. A 















TABLE IV. - COMPARISON BETWEEN MEASURED AND PREDICTED EFFICIENCIES AT DESIGN 

TEST PRESSURE RATIO 



Turbine 

76 

51 

42 A 

42 B 

25 U 

25 T 

12 

13-Cyl 

13-Cone 

% equivalent speed 

100 

100 

100 

100 

100 

100 

100 

100 

100 

Measured efficiency 

.922 

.895 

.920 

.911 

.874 

.882 

.863 

.845 

.855 

Predicted efficiency 

.923 

.907 

.916 

.915 

.881 

.876 

.865 

.850 

.853 

% equivalent speed 

90 

110 

no 

no 

no 

no 

no 

no 

no 

Measured efficiency 

.904 

.901 

.921 

.912 

.884 

.892 

.876 

.857 

.863 

Predicted efficiency 

.900 

.923 

.919 

.910 

.878 

.890 

.877 

.861 

.865 

% equivalent speed 

80 

80 

90 

90 

80 

80 

80 

80 

80 

Measured efficiency 

.884 

.845 

.906 

.901 

.836 

.860 

.821 

.799 

.808 

Predicted efficiency 

.869 

.848 

.902 

.913 

.840 

.845 

.825 

.795 

.804 















TABLE V. - NEW DEFINITIONS OF INPUT VARIABLES 
FOR ELLIPTICAL LEADING AND TRAILING EDGES 


RI1 Semi-major axis of leading-edge 
ellipse 

RI2 Semi-minor axis of leading-edge 
el 1 ipse 

BETI1 

Angle of major axis of leading-edge 
el 1 ipse. 


BETI2 

Not used 

MSP 1(1) 

m-coordinate of tangent point of 

spline curve on leading-edge 
ellipse for upper surface 
(Surface 1). Previously not 
required. 


MSP2(1) 

Same as MSP1(1), but for lower 
surface (Surface 2). 

R01 Semi-major axis of trailing edge 
el 1 ipse 

R02 Semi-minor axis of trailing edge 
el 1 ipse 

BET01 

Angle of major axis of trailing 
edge ellipse 


BET02 
Not used 

MSPl (SPLN01 ) 

Same as MSPi(l), but at trailing 
edge 

MSP2(SPLN02) 

Same as MSPl (SPLN01) , but for lower 
surface 

Note: If the trailing edge is circular, but the 

leading edge is elliptical, it is necessary to 
specify that the values of MSPl(SPLNOl) and 
MSP2(SPLN02) . 



START 



Figure 1. - Diagram of procedure used to calculate stage efficiency. 













k SOD S A l\ 


SECONDARY LOSS CORRELATIONS 

MORRIS AND HOARE 
DUNHAM CORRELATION 


EXPERIMENTAL DATA GIVEN IN REF. 19 

□ MEASUREMENTS INCLUDED DOWNSTREAM WALL 
BOUNDARY LAYER 

A MEASUREMENTS MAY HAVE INCLUDED 
DOWNSTREAM WALL BOUNDARY LAYER 
O MEASUREMENTS EXCLUDED DOWNSTREAM 
WALL BOUNDARY LAYER 


a 4 

.03 

oo 

O 


l M 

.02 

i « 

.01 


19 POINTS 


-TOTAL SECONDARY LOSS 

> SECONDARY LOSS DUE TO 
INLET BOUNDARY LAYER 

20 POINTS~ EQLJATION 2 

J l 



0 .01 .02 .03 .04 .05 

UPSTREAM WALL BOUNDARY LAYER 
DISPLACEMENT THICKNESS/CHORD. 


6^/c 


Figure 2. - Comparison of secondary loss correlations with 
the experimental data given in reference 19. 



RATIO 

Figure 3. - Comparison of measured and 
predicted efficiencies for stage 76. 
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Figure 4. - Comparison of measured and predicted 
stage efficiency for stage 51. 
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Figure 5. - Comparison of measured and predicted 
efficiencies for stage 42 A. 
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Figure 6. - Comparison of measured and predicted 
efficiencies for stage 42B. 
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Figure 7. - Comparison of measured and predicted 
stage efficiency for stage 25 U. 
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Figure 8. - Comparison of measured and predicted 
stage efficiency for stage 25 T. 
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Figure 9. - Comparison of measured and predicted ' 
stage efficiency for stage 12. 



Figure 10. - Comparison between measured and pre- 
dicted stage efficiency for turbine stage 13 cyl. 
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Figure 11. - Comparison between measured and predicted 
stage efficiency for turbine stage 13-Cone. 
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(b) Constant stage pressure ratio. 

Figure 12. - Comparison between measured 
and predicted stage efficiency at different 
clearances for stage 13-cyl. 
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Figure 13. - Loss distribution for each stage at design condition. 


PERCENT OF TOTAL LOSS 











ROTOR EXIT FLOW ANGLE, deg 




Figure 14. - Comparison of measured and 
predicted rotor exit flow angles for three 
turbines at the design point. 


Figure 15. - Schematic of procedure used for trailing edge ejection. 











Figure 16. - Schematic of nomenclature used for elliptical 
blade leading edge. 
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Figure 17. - Comparison of grids generated starting at different 
spanwise locations. 



Figure 18. - Sketch of path used to determine stage 
efficiency. 
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